The future of biomolecular simulation in the pharmaceutical industry: what we can learn from aerodynamics modelling and weather prediction. Part 1. understanding the physical and computational complexity of in silico drug design

The predictive power of simulation has become embedded in the infrastructure of modern economies. Computer-aided design is ubiquitous throughout industry. In aeronautical engineering, built infrastructure and materials manufacturing, simulations are routinely used to compute the performance of potential designs before construction. The ability to predict the behaviour of products is a driver of innovation by reducing the cost barrier to new designs, but also because radically novel ideas can be piloted with relatively little risk. Accurate weather forecasting is essential to guide domestic and military flight paths, and therefore the underpinning simulations are critical enough to have implications for national security. However, in the pharmaceutical and biotechnological industries, the application of computer simulations remains limited by the capabilities of the technology with respect to the complexity of molecular biology and human physiology. Over the last 30 years, molecular-modelling tools have gradually gained a degree of acceptance in the pharmaceutical industry. Drug discovery has begun to benefit from physics-based simulations. While such simulations have great potential for improved molecular design, much scepticism remains about their value. The motivations for such reservations in industry and areas where simulations show promise for efficiency gains in preclinical research are discussed. In this, the first of two complementary papers, the scientific and technical progress that needs to be made to improve the predictive power of biomolecular simulations, and how this might be achieved, is firstly discussed (Part 1). In Part 2, the status of computer simulations in pharma is contrasted with aerodynamics modelling and weather forecasting, and comments are made on the cultural changes needed for equivalent computational technologies to become integrated into life-science industries

Modelling the binding mode of macrocycles: Docking and conformational sampling

Drug discovery is increasingly tackling challenging protein binding sites regarding molecular recognition and druggability, including shallow and solvent-exposed protein-protein interaction interfaces. Macrocycles are emerging as promising chemotypes to modulate such sites. Despite their chemical complexity, macrocycles comprise important drugs and offer advantages compared to non-cyclic analogs, hence the recent impetus in the medicinal chemistry of macrocycles. Elaboration of macrocycles, or constituent fragments, can strongly benefit from knowledge of their binding mode to a target. When such information from X-ray crystallography is elusive, computational docking can provide working models. However, few studies have explored docking protocols for macrocycles, since conventional docking methods struggle with the conformational complexity of macrocycles, and also potentially with the shallower topology of their binding sites. Indeed, macrocycle binding mode prediction with the mainstream docking software GOLD has hardly been explored. Here, we present an in-depth study of macrocycle docking with GOLD and the ChemPLP scores. First, we summarize the thorough curation of a test set of 41 protein-macrocycle X-ray structures, raising the issue of lattice contacts with such systems. Rigid docking of the known bioactive conformers was successful (three top ranked poses) for 92.7% of the systems, in absence of crystallographic waters. Thus, without conformational search issues, scoring performed well. However, docking success dropped to 29.3% with the GOLD built-in conformational search. Yet, the success rate doubled to 58.5% when GOLD was supplied with extensive conformer ensembles docked rigidly. The reasons for failure, sampling or scoring, were analyzed, exemplified with particular cases. Overall, binding mode prediction of macrocycles remains challenging, but can be much improved with tailored protocols. The analysis of the interplay between conformational sampling and docking will be relevant to the prospective modelling of macrocycles in general

Analysis of Arm Movement Strategy in Virtual Catching Task

In this paper, we explored how the arm movement pattern as well as the related strategy of the children with Cerebral Palsy (CP) and the healthy children can be changed in the virtual catching task on a previously proposed rehabilitation system. We recruited 50 healthy children from elementary school, and 3 children with CP as subjects to classify their arm movement pattern/strategy. As a result of the classification, we identified three arm movement stages : Initial position, Reaching path, and Waving form, as well as movement pattern strategy under each movement stage. Based on the classified pattern, we compared the differences in the time series changes of movement strategy between healthy children and the children with CP. The results show there is a significant difference in the strategy of arm movements in the Initial position between healthy and CP children

The future of biomolecular simulation in the pharmaceutical industry: what we can learn from aerodynamics modelling and weather prediction. Part 1. understanding the physical and computational complexity of in silico drug design

The predictive power of simulation has become embedded in the infrastructure of modern economies. Computer-aided design is ubiquitous throughout industry. In aeronautical engineering, built infrastructure and materials manufacturing, simulations are routinely used to compute the performance of potential designs before construction. The ability to predict the behaviour of products is a driver of innovation by reducing the cost barrier to new designs, but also because radically novel ideas can be piloted with relatively little risk. Accurate weather forecasting is essential to guide domestic and military flight paths, and therefore the underpinning simulations are critical enough to have implications for national security. However, in the pharmaceutical and biotechnological industries, the application of computer simulations remains limited by the capabilities of the technology with respect to the complexity of molecular biology and human physiology. Over the last 30 years, molecular-modelling tools have gradually gained a degree of acceptance in the pharmaceutical industry. Drug discovery has begun to benefit from physics-based simulations. While such simulations have great potential for improved molecular design, much scepticism remains about their value. The motivations for such reservations in industry and areas where simulations show promise for efficiency gains in preclinical research are discussed. In this, the first of two complementary papers, the scientific and technical progress that needs to be made to improve the predictive power of biomolecular simulations, and how this might be achieved, is firstly discussed (Part 1). In Part 2, the status of computer simulations in pharma is contrasted with aerodynamics modelling and weather forecasting, and comments are made on the cultural changes needed for equivalent computational technologies to become integrated into life-science industries

Exploring the conformational dynamics of alanine dipeptide in solution subjected to an external electric field: A nonequilibrium molecular dynamics simulation

In this paper, we investigate the conformational dynamics of alanine dipeptide under an external electric field by nonequilibrium molecular dynamics simulation. We consider the case of a constant and of an oscillatory field. In this context we propose a procedure to implement the temperature control, which removes the irrelevant thermal effects of the field. For the constant field different time-scales are identified in the conformational, dipole moment, and orientational dynamics. Moreover, we prove that the solvent structure only marginally changes when the external field is switched on. In the case of oscillatory field, the conformational changes are shown to be as strong as in the previous case, and non-trivial nonequilibrium circular paths in the conformation space are revealed by calculating the integrated net probability fluxes.Comment: 23 pages, 12 figure

Facile heterocyclic synthesis and antimicrobial activity of polysubstituted and condensed pyrazolopyranopyrimidine and pyrazolopyranotriazine derivatives

Reaction of 6-amino-3-methyl-4-(substituted phenyl)-1,4-dihydropyrano[2,3-c]pyrazole-5-carbonitrile (1) with triethylorthoformate followed by treatment with hydrazine hydrate, formic acid, acetic acid, phenylisocyanate, ammonium thiocyanate and formamide afforded the corresponding pyranopyrimidine derivatives 2–6. Cyclocondensation of 1 with cyclohexanone afforded pyrazolopyranoquinoline 7. One-pot process of diazotation and de-diazochlorination of 1 afforded pyrazolopyranotriazine derivative 8, which upon treatment with secondary amines afforded 9 and 10a-c. Condensation of 2 with aromatic aldehyde gave the corresponding Schiff bases 11a,b, the oxidative cyclization of the hydrazone with appropriate oxidant afforded 11-(4-fluorophenyl))-2-(4-substitutedphenyl)-10-methyl-8,11-dihydropyrazolo-[4\u27,3\u27:5,6]pyrano[3,2-e][1,2,4]triazolo[1,5-c]pyrimidines (12a,b). Structures of the synthesized compounds were confirmed by spectral data and elemental analysis. All synthesized compounds were evaluated for antibacterial and antifungal activities compared to norfloxacin and fluconazole as standard drugs. Compounds 9, 10c, 12a and 15 were found to be the most potent antibacterial agents, with activity equal to that of norfloxacin. On the other hand, compound 5 exhibited higher antifungal activity compared to fluconazole

How Thioredoxin Dissociates Its Mixed Disulfide

The dissociation mechanism of the thioredoxin (Trx) mixed disulfide complexes is unknown and has been debated for more than twenty years. Specifically, opposing arguments for the activation of the nucleophilic cysteine as a thiolate during the dissociation of the complex have been put forward. As a key model, the complex between Trx and its endogenous substrate, arsenate reductase (ArsC), was used. In this structure, a Cys29Trx-Cys89ArsC intermediate disulfide is formed by the nucleophilic attack of Cys29Trx on the exposed Cys82ArsC-Cys89ArsC in oxidized ArsC. With theoretical reactivity analysis, molecular dynamics simulations, and biochemical complex formation experiments with Cys-mutants, Trx mixed disulfide dissociation was studied. We observed that the conformational changes around the intermediate disulfide bring Cys32Trx in contact with Cys29Trx. Cys32Trx is activated for its nucleophilic attack by hydrogen bonds, and Cys32Trx is found to be more reactive than Cys82ArsC. Additionally, Cys32Trx directs its nucleophilic attack on the more susceptible Cys29Trx and not on Cys89ArsC. This multidisciplinary approach provides fresh insights into a universal thiol/disulfide exchange reaction mechanism that results in reduced substrate and oxidized Trx

Excitons in a Photosynthetic Light-Harvesting System: A Combined Molecular Dynamics/Quantum Chemistry and Polaron Model Study

The dynamics of pigment-pigment and pigment-protein interactions in light-harvesting complexes is studied with a novel approach which combines molecular dynamics (MD) simulations with quantum chemistry (QC) calculations. The MD simulations of an LH-II complex, solvated and embedded in a lipid bilayer at physiological conditions (with total system size of 87,055 atoms) revealed a pathway of a water molecule into the B800 binding site, as well as increased dimerization within the B850 BChl ring, as compared to the dimerization found for the crystal structure. The fluctuations of pigment (B850 BChl) excitation energies, as a function of time, were determined via ab initio QC calculations based on the geometries that emerged from the MD simulations. From the results of these calculations we constructed a time-dependent Hamiltonian of the B850 exciton system from which we determined the linear absorption spectrum. Finally, a polaron model is introduced to describe quantum mechanically both the excitonic and vibrational (phonon) degrees of freedom. The exciton-phonon coupling that enters into the polaron model, and the corresponding phonon spectral function are derived from the MD/QC simulations. It is demonstrated that, in the framework of the polaron model, the absorption spectrum of the B850 excitons can be calculated from the autocorrelation function of the excitation energies of individual BChls, which is readily available from the combined MD/QC simulations. The obtained result is in good agreement with the experimentally measured absorption spectrum.Comment: REVTeX3.1, 23 pages, 13 (EPS) figures included. A high quality PDF file of the paper is available at http://www.ks.uiuc.edu/Publications/Papers/PDF/DAMJ2001/DAMJ2001.pd

Multiple pH Regime Molecular Dynamics Simulation for pK Calculations

Ionisation equilibria in proteins are influenced by conformational flexibility, which can in principle be accounted for by molecular dynamics simulation. One problem in this method is the bias arising from the fixed protonation state during the simulation. Its effect is mostly exhibited when the ionisation behaviour of the titratable groups is extrapolated to pH regions where the predetermined protonation state of the protein may not be statistically relevant, leading to conformational sampling that is not representative of the true state. In this work we consider a simple approach which can essentially reduce this problem. Three molecular dynamics structure sets are generated, each with a different protonation state of the protein molecule expected to be relevant at three pH regions, and pK calculations from the three sets are combined to predict pK over the entire pH range of interest. This multiple pH molecular dynamics approach was tested on the GCN4 leucine zipper, a protein for which a full data set of experimental data is available. The pK values were predicted with a mean deviation from the experimental data of 0.29 pH units, and with a precision of 0.13 pH units, evaluated on the basis of equivalent sites in the dimeric GCN4 leucine zipper